Light mass galileons: Cosmological dynamics, mass screening and observational 

constraints 
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In this paper, we examine the cosmological viability of a light mass galileon field consistent 
with local gravity constraints. The minimal, L3 = □0(9 M 0) 2 , massless galileon field requires an 
additional term in order to give rise to a viable ghost free late time acceleration of Universe. The 
desired cosmological dynamics can either be achieved by incorporating an additional terms in the 
action such as (£4, L5) — the higher order galileon Lagrangians or by considering a light mass 
field a la galileon field potential. We analyse the second possibility and find that: (I) The model 
produces a viable cosmology in the regime where the non- linear galileon field is subdominant, (2) 
The Vainshtein mechanism operates at small scales where the non-linear effects become important 
and contribution of the field potential ceases to be significant. Also the small mass of the field under 
consideration is protected against strong quantum corrections thereby providing quantum stability 
to the system. 
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I. INTRODUCTION 

Modified theories of gravity pose a serious alternative 
to dark energy, an exotic cosmic fluid, needed to account 
for late time cosmic acceleration in the framework of 
standard lore[l|-0] ■ The task of building an alternative to 
Einstein's theory is very challenging as the latter fits with 
the observation with a great accuracy locally thereby a 
large scale modification is either felt locally or the frame- 
work gets reduced to AC DM. A viable alternative theory 
of gravity should satisfy important requirements: (1) It 
should be close to ACDM but yet distinguishable from 
it, (2) Theory should be free from ghost and tachyon in- 
stabilities and (3) Theory should not conflict with local 
physics. The third requirement is often very stringent 
and its compliance needs the invoking of special mech- 
anisms. The modified theories necessarily include addi- 
tional scalar degree(s) of freedom in some form or the 
other. And to be relevant to late time cosmic accelera- 
tion, it should be light mass entity which, on the other 
hand, could cause a havoc as a fifth force never seen in 
the laboratory or in the solar neighborhood. It therefore 
necessary to screen out the effect of the fifth force in a 
delicate manner. 

Broadly, there are two methods of mass screening, the 
chameleon mechanism and the Vainshtein effect. The 
chameleon scenario 0] relies on the direct coupling of 
matter to scalar field with potential. The mass of the 
field becomes dependent on the density of environment 
which justifies the designation "chameleon" for such a 
field. The chameleon potential is chosen such that the 
effective mass of the field increases with local matter den- 
sity thereby leading to suppression of fifth force. The 
chameleon mechanism is an extremely powerful tool for 
mass screening but as noticed by many authors (see e.g. 



[Io| ) the numerical integration of the system for a spher- 
ically symmetric background hits singularity because of 
the form of the chameleon potential. Hence an extreme 
fine-tuning of the initial conditions is required to avoid 
the problem. The chameleon theories arc also plagued 
with the problem of large quantum corrections because 
of the large mass of the chameleon field required to pass 
the local gravity tests [ll| . 

The Vainshtein mechanism [l2| is a superior field theo- 
retic method of mass screening. It was invented by Vain- 
shtein in 1972 to address the discontinuity problem in 
massive gravity of Pauli-Ficrz. It relics on non-linear 
derivative term of the type £3 = (d^ip^Dcj) where (j> is 
scalar degree of freedom of hclicity zero graviton in this 
case. The dynamics of non-linear term gives rise to a 
miraculous phenomenon: Around a massive body, in a 
large radius dubbed Vainshtein radius, fifth force is sup- 
pressed switching off any modification to gravity locally. 
The field is strongly coupled to itself and hence becomes 
weakly coupled to matter sector. The DGP model [l3[ 
contains such a non-linear term in the so called decou- 
pling limit responsible for the compliance of the model 
with local gravity constraints The non-standard 

kinetic term also occurs in Kaluza-Klcin reduction of 
Gauss-Bonnet gravity to four dimensional space time (l5j 
which makes clear why the theory is free from Ostrograd- 
sky ghosts. The role of the scalar degree of freedom is 
played by dilaton in this case. The field <f> due to the 
presence of an underlying symmetry in flat space time 
was termed as galileon. There exist higher order galileon 
Lagrangians L4 and L5 that contain higher order non lin- 
ear derivative terms in which four and five galileon fields 
participate respectively. Recently, more general galileon 
action were constructed 1 171 an d their cosmological 
implications were investigated [18[. They belong to a 
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more general class of models first introduced in [19( . 

The galileon field with the lower order term L3 is suffi- 
cient to take care of the local gravity constraints whereas 
L5 does not contribute to mass screening and though L4 
can effect the numerics of Vainstcin radius but adds noth- 
ing to the underlying physics of Vainshtcin mechanism. 

However, galileon system with L3 term alone can not 
give rise to late time acceleration of universe. It was 
first demonstrated i n 1201 that at least one higher order 
Lagrangian, say L4 |23( be added to action in order to 
produce a stable de Sitter solution. This is the analogue 
of the DGP model, where the self-accelerating branch is 
unstable because of the presence of ghost 24[ . Also phys- 
ical implications of the cubic Galileon term coupled non- 
minimally to the metric was first studied in [2ll |22| . Since 
we are working in a phenomenological setting, we could 
restrict ourselves to L3 but add a potential to galileon to 
produce late time cosmic acceleration. We do not assign 
a field theoretic mechanism to produce the potential for 
galileon, perhaps one would need some non-perturbative 
machinery to do the job. 

Our proposal can be seen as an attempt to re- 
place chameleon mechanism by the Vainshtein effect: 
In chameleon scenario, the choice of specific form of 
chameleon potential is crucial to the model which needs 
to be extremely fine-tuned locally for the successful im- 
plementation of the underlying chameleon ideology. 



II. BASIC SETUP AND EQUATIONS OF 
MOTION 



It is well known that in General theory of Relativity, 
the Newtonian and longitudinal gravitational potentials 
are equal provided that anisotropic stress tensor is ab- 
sent. Modified gravity theories usually predict a differ- 
ence between these two potentials — a gravitational slip. 
At the linear level of perturbations, the modified theories 
of gravity give rise to an anisotropic stress tensor which 
is really not a signature of a particular modified theory 
as usually claimed. At the linear level of perturbations, 
the modified theories of gravity give rise to an anisotropic 
stress tensor which is really not a distinctive signature of 
modification of gravity as usually claimed. A modified 
theory of gravity can, in general, be described in terms 
of an imperfect fluid giving rise to a gravitational slip. 
Thus we can not distinguish whether the slip is caused 
by the presence of an imperfect fluid or by the modifica- 
tion of gravity. Since the observations, at present, do not 
show any difference between the gravitational potentials 
(see e.g. HH), it is interesting to study models that do 
not produce a gravitational slip. 

Following the results derived in [26|, we deduce that 
the most general action with second order differential 
equation for the metric and the scalar field which do not 
produce a slip is of one of the following form 
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F((j))R+ ( F' (<f>) In X + G(<j>))a(j) 
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or 



So = 




G{4>,X)D<f> + K(<f>,X) 



XR-G 



^^)] (2) 



where X 



The two actions are different but in both cases we do 
not have a gravitational slip. The second action S2 is a 
generalization of the KGB model studied in [27j . We also 
note that for S\ in the Brans-Dicke form <f>R, the presence 
of \nX Dtp is important to cancel the slip between the 
two potentials. 

Keeping in mind a well motivated model which appears 
in the decoupling limit of DGP, we shall consider the 
second action hereafter. We also introduce a coupling of 
field to matter, therefore we work in the Einstein frame 
where the theory does not gives rise to a slip. Obviously 
in the Jordan frame the same theory would produce a slip 
because of the conformal factor as noticed in J28[. The 
conformal transformation would produce terms which are 
not in the form of the original action. 

The second action is the most general action free from 
Ostrogradsky ghost problem, it does not produce a grav- 
itational slip and contains the decoupling limit of DGP 
as a sub-class. But the first action can be seen as a mod- 
ified Brans-Dicke model which can produce interesting 
local and cosmological solutions. 

In what follows, we shall consider a simple model with 
G = X which corresponds to the decoupling limit of 
DGP, F = and K = X - V{4>). With these choices, 
our action has the form (in Einstein frame) 
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(3) 



where = 8irG is the reduced Planck mass, M is a 
energy scale, (a, /?) are dimensionless constants and V is 
the potential for the field. 

This action corresponds to the coupled quintessence field 
with a galileon like correction (V0) 2 D<j>. Variation of the 
action ^ gives the following equations of motion 
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where 



a 
M 3 



(6) 



where ' denotes the derivative wrt </>. 

In the discussion to follow, we shall investigate the 
cosmological dynamics and local gravity constraints in 
the model based upon the equations of motion §5$ . 



III. FIFTH FORCE, SCREENING EFFECT AND 
LOCAL PHYSICS 



the potential is dominant and this part would be stud- 
ied in the next section. We note there is also an in- 
termediate regime where the potential is of the order of 
L 3 . We would not consider the analysis of this regime 
here. In fact at very large distances, the approximations 
made in this section can not be trusted, the expansion 
of the Universe should be taken into account too. Last 
but not least, as there is no physical system isolated at 
large scales, it is not known how the screening mecha- 
nism would operate for N particles if each one of them 
possess its own Vainshtein radius — the so-called elephant 
problem. 

For the scales of interest mentioned above, the fifth 
force is 



In our model, in general, we might have both the 
screening mechanisms operating due to the presence of 
potential and the galileon term. We wish to arrange 
things in such a way that potential contributes negligi- 
bly to local physics giving way to Vainshtein mechanism. 
The potential should come into action at large scales giv- 
ing rise to cosmic acceleration. 

We now consider spherically symmetric, static solution 
of equation (J5J) . In flat Minkowski background and in 
presence of non relativistic matter, we have 



-A [rV(r)] 
H dr 



2a d 
M 3 r 2 dr 
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M, 



-Pm 

(7) 
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We consider that the same phenomenon is at the ori- 
gin of late time acceleration and local suppression of fifth 
force in the action. We therefore associate to both the 
non-linear term and the potential the same energy scale 
M. We have <j> ~ M and V — M 4 . As we previously men- 
tioned, the chameleon mechanism generates huge fine- 
tuning conditions. Thus if we want to prefer the Vain- 
shtein mechanism, we should demand that at small scales 
3P7?3F V<P'{rf] » ™. Then considering a ~ 1 we 
may have a negligible potential contribution for scales 
L <C 1/M. This scale should at least be larger that the 
apogee which tells us that M <C 10~ 15 cV. Also we should 
assume that the potential term is negligible compared to 
the matter part. Considering p m ~ lg/cm~ 3 , we have 
M <C 10~ 3 eV. This condition could always be verified 
as soon as the first one is fixed. We shall see later that 
for viable cosmology, the potential should be of the order 
of critical density which implies M ~ 10 _3 cV. It seems 
contradictory with the previous statement. In fact, from 
the energy scale M, the results will also depend on the 
form of the potential. We found that for inverse power 
law potentials, we are in the chameleon regime and hence 
divergences can occur. We will take up this issue in the 
next section where we shall consider the class of poten- 
tials in detail which give negligible contribution locally 
for M ~ l(T 3 cV 

In that case, neglecting the potential, the system re- 
duces to the standard Galileon with L3. At large scales, 



/3M 3 r 
AaM p i 



r / 



(8) 



where ry is the Vainshtein radius of the body with 
mass M and r s is the Schwarzschild radius. 



8afi M pl r s 
M 3 



(9) 



The Vainshtein radius is defined as the distance from 
where onwards the fifth force becomes comparable with 
the gravitational force. Considering the Sun, we need 
a Vainshtein radius larger than the size of solar system 
which implies that M < 1(T 7 (a/3) 1/3 cV. The scale M 
plays a fundamental role in local physics determined by 
Vainshtein mechanism. This is also the scale associated 
with the potential thereby to the cosmology of the model. 
In order to get late time acceleration of the Universe, we 
need to fix, M ~ 10 -3 eV. On the other hand, in order to 
have Vainshtein radius larger than the solar system, we 
should demand that aj3 > 10 12 . Another option would be 
to choose a ~ 1 or M 4 <C A where A is the cosmological 
constant. This would make the model even less natural 
than the cosmological constant in terms of an unnatural 
value. Perhaps, it points towards a possibility that the 
scale could be linked to some large extra dimension a la 
DGP model. Finally, one could also consider two dif- 
ferent scales (Mi,M2) in the model associated with the 
galileon term and the potential respectively. We shall not 
consider this possibility here. 

In our analysis we have considered two different poten- 
tials 



V(<t>) = M 4 



e M P i 



_ V2 4> 

g M pl 



and 



V(<j>) = M 4 



cosh 



7^ 



(10) 



(11) 
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where A, p, 7 are constants and m is a positive integer. 
In both cases, we found numerically that the potential 
is negligible locally and that the eq.© can be approxi- 
mated by a vanishing potential. 

The system appears as the decoupling limit of DGP and 
as it is well known that the Vainshtein mechanisms pro- 
tects the model in a static spherically symmetric config- 
uration. In the next section we shall demonstrate that 
the potential is dominant at large scales and gives rise to 
a correct cosmology contrary to the same model without 
potential [Io|. 



IV. ON THE BOUND FROM THE TIME 
VARIATION OF G 

All the theories described in the Einstein frame by a 
time-varying coupling to matter leads to a variation of 
the effective Newton constant (G) as seen in the Jor- 
dan frame. There are two effects which modify G p9| : 
The exchange of hclicty-0 modes which is suppressed 
in our model because of the Vainshtein mechanism and 
the rescaling of the coordinates because of the conformal 
transformation. In our case the time variation of G oc- 
curs because of the rescaling factor only. Hence it is easy 
to see that 



\G/G\ m /3<j>/M pl . 



(12) 



It was shown in Ref.[23| that a model involving only 
derivatives of the fields would give 



\G/G\ today 



(13) 



But the observational constraints from the Lunar Laser 
Ranging (LLR) gives |G/G| to day < 0.01i? (H which 
translates into a restriction on /3, namely, j3 < 0.01. In 
the presence of a potential, the situation is drastically 
modified. In fact, in this case, we can assume that the 
field sits at the minimum of the effective potential. As a 
result, in case of the potential (flU)) . we have 



\G/G\ 



today 



6/3 2 



Ml 



— 2^0 

P2 



Hence the LLR constraint gives 



< 



vJ4 



p\ 



25 



(14) 



(15) 



In the next section, we will consider cosmology for (pi = 
20 /12) which means that (3 < 1. Therefore we see that 
the LLR bound will not bring additional constraints as 
soon as the standard cosmological constraints are satis- 
fied. 



V. COSMOLOGICAL DYNAMICS OF THE 
MODEL 

In a spatially flat FLRW background, the equations of 
motion take the form 



12 /? 

3M 2 ,H 2 = Pm + Pr + t-(l- JL H( j> j + V{4>) 

(16) 



Af 2 l( 2F + 3ff 2 ) = -^-^(l + ^>) 



3a • 



(17) 



M p r" 1 - ■ — - M 3 



-^-p m =<t> + 3H<f>- —cf,( m z <t> + #</> + 2H<t> 

(18) 



The equation of conservation which can be derived 
from the previous equations is 

P i 



Pm + 3i? Pm — 



(19) 

Pr + AHp r = (20) 

Let us introduce the following dimensionless quantities 



(21) 



(22) 



needed to cast the evolution equations in the form of 
an autonomous system 



d.T _ / (j) H \ 

dy l /3. 

dN 



de _ / H\ 
dN _e V^ + Jp) 

dO r / H 

^ = VfaA'(i-r) 



(23) 

(24) 

(25) 

(26) 
(27) 



where TV = In a, T 



and 



H _ 2(1 + e)(-3 + 3y 2 - fl r ) - 3x 2 (2 + 4e + e 2 ) 
H 2 ~ 4 + 4e + £ 2 e 2 

, V6xe(y 2 \- pn m ) 



4 + 4e + x 2 e 2 



(28) 



3a; e — x 



(l2 + e(3 + 3y 2 - a )) + 2\/% 2 A - pn m ) 



H<t> 



x(A + 4e + x 2 e 2 ) 



(29) 



In the case of an exponential form of the potential, 
F = 1 and A is a constant. Therefore the system reduces 
to the set of four equations as usual. 
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For e = 0, we recover the standard coupled dark-energy 
model |3l| . The system has the same dynamical phase 
plane as the coupled quintessence except one additional 



de-Sitter solution, for (e = —2, A = 0, y 



1). This 



solution exists only if the equation A 2 T(A) = admits the 
solution A = 0. Before getting to numerical analysis, we 
should broadly summarize the constraints on the model 
parameters, we have parameters (/zi, ^2) in the potential 
(jTUJ) which control the different phases of the dynamics. 
One of the parameters /Xi , gives rise to the radiation era 
whereas, (12 is responsible for the late time cosmic accel- 
eration. The matter era is principally controlled by /?. 
In order to have viable thermal history (Q, r = 1 — 4//zf), 
we need to impose a constraint on fix. Indeed, the nucle- 
osynthesis demands that f2^ = 4/yU 2 < 0.01 during the 
radiation era which imposes restriction on fi± , namely, 
\X\ > 20. For stable acceleration of the Universe we need 
to impose constraints on other parameters such that, 
/i2 < v2 and also (3 < (3 — /zf)//^. During the Dark 
Energy phase, uj c s = — 1 + /z§/3. We shall make use of 
this information while carrying out the numerical analy- 
sis. Indeed, in the numerics to follow, we will take small 
values for (3 and adhere to (zti,zt 2 ) = (20,0.5). During 
the evolution, A would change from \i\ to \i2 hence if \i2 is 
small enough, it would imply a small variation of the po- 
tential which would give rise to quasi-dc-Sittcr universe, 
w off ~ -1. 

We have numerically analysed the autonomous system 
and find that the cosmology of the model is similar to 
the case of coupled quintessence Fig. (|l|2|) . The influence 
of the higher derivative term disappears fast with evolu- 
tion leaving behind the coupled quintessence. We have 
not listed the various fixed points and their correspond- 
ing eigenvalues to avoid the repetition of the earlier work 
on the similar theme. We have plotted the dimension- 
less density parameter and field energy density versus 
the scale factor on the log scale for the double exponen- 
tial potential. The plots show in case of the overshoot a 
tracker behaviour. The parameters can be conveniently 
fixed to produce a viable late time cosmology. 



VI. COSMOLOGICAL PERTURBATION 

In this section we analyse the perturbations of the 
model. As mentioned in the introduction, the action is 
built in such away that the two gravitational potentials 
are same. Hence we consider the following metric 

ds 2 = -(1 + 2^)dt 2 + a(t) 2 (l - 2^)dx 2 (30) 
In the subhorizon approximation, we have 



S m + (2H 



J_ 

Mr,, 



^GeffPmSm = (31) 
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FIG. 1: (Top): Density parameters f2 of matter, radiation and 
galileon field for potential (|10|) are shown here with /ii = 20 

, fJ.2 = 0.5 and /? = 0.1. (Bottom): Density log ^ , iM § [ri ^ for 

the 3 fields and the same parameters (/zi, /Z2, /3). 
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FIG. 2: The cosmic evolution of the field equation of state 
w</> and the total effective equation of state to e ff for fj,i = 20 , 
pa = 0.5 and /3 = 0.1. 



Where, 



G 



eff 



G 



f_s_ 

\M 3 M pl 



4> 2 + 2f3 



2 - 



M 6 M' A , 



(32) 

where the gauge-invariant density contrast is defined 
as 5 m — + 3Hv and v is the velocity of the fluid. 

In Fig.©, we show the evolution of G e ff as a func- 
tion of the redshift. Because of the non-linear term of 
the galileon field, we have at large redshift G e g similar 
to General Relativity which is consistent with BBN con- 
straints, whereas G c g = G(l + 2/3 2 ) > G in the cou- 
pled quintessence model. Thus the non-linear term of 
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FIG. 3: Evolution of GW/G for the potential (fT0|) with fii = 
20, fi 2 = 0.5 and = 0.5. 

the galileon field can alleviate constraints on the model. 
We have seen that the bound from the LLR is weak, 
hence we should expect to have large values of the cou- 
pling /3. However, we shall notice in the section to follow 
that the model under consideration reduces to coupled 
quintessence scenario at large scales thereby giving rise 
to strong cosmological constraints on /3 as first noticed 
in Ref.(3l|. 

VII. DATA ANALYSIS 
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FIG. 4: Evolution of the growth factor / for potential ({TO}) 
and respectively. The integration is done from Q m ~ 0.98 
for both the cases with f3 = 0.1. For the potential (|10|l we 
have taken fii — 20 and /12 = 0.5 and for the potential (|11[1 
we have taken 7 = 30 and m = 0.5. The red line corresponds 
to ACDM model. 



To get some insight into the parameter window for the 
models, we constrain them using Supernovae, BAO and 
CMB data. We have also used the data of the growth 
factor / [32[ defined as 



/ 



d In 8 m 
din a 



(33) 



We show the evolution of the growth factor / for the 
two potentials studied in this paper. In both cases, 
the deviation from ACDM model is not significant. Wc 
should also mention that / > 1 at z = 3 which means 
that the growth is faster than in Einstein-de-Sitter model. 
This is opposite to the prediction of Dark Energy sce- 
nario. This sounds very strange but since the deviation 
from ACDM is not significant, the latter does not effect 
the statistical analysis. Only future observations can con- 
form this result. 

We define 



For the SN type la observation we consider Union2.1 
data compilation [33| of 580 data points. The observable 
for the SN type la is the distance modulus \i which is 
defined as, \i = m — M = 5 log Dl + /j,q, where m and M 
are the apparent and absolute magnitudes of the Super- 
novae and luminosity distance Dl,{z) = (1 + z) J Q Z 

and jUo = 5 log ^ +25 is a nuisance parameter and 

should be marginalized. For this observation the \ 2 is 
defined as, 



XsnWo, Q) = 2-^ Z t„\2 ( 36 ) 

t=l 



2 2 2 2 2 / \ 

X = XGrowth + XsN + XBAO + XcMB- (34) 



where y 2 - 



rowth 



is defined as, 



Where, is the uncertainty in the distance modulus. 
After marginalizing /iq as [3~i| we can have, 



Xgt 



uth 



(») = £ 



fobs(Zj) - fth{zi,t 
CTf(Zi) 



(35) 



Where, / b s is the observational data and o $ is the la 
error. f t h is the theoretically calculated value of growth 
factor from the eqn (|33p . 8 is the model parameter. 



X 2 sn(0)=A 



c 



(37) 



Where, 
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Xbao 1S defined as, 



i=l 
580 



b{9) = y: 



0>(Zj) 2 



= 1 

580 



(38) 
(39) 
(40) 



For BAO we have used the data of rwi 2 ^ ^ [H, 
where z+ is the decoupling time ~ 1091, the co- 
moving angular-diameter distance (Ia(z) = Jq 2 "^pj an d 



Dv(z) = \ cIa(z) jjj^j) ■ Required data for this analysis 
is given in table [7] 



Where, 



X 



BAO 



D v (0.106) 

ZV(0.2) 
<Ia{z*) 
D v (0.35) 

IV (0.44) 

Dv 0.6) 

dA{Zt,) 

IV (0.73) 



- 30.95 \ 

- 17.55 

- 10.11 
-8.44 

- 6.69 

- 5.45 



(41) 



(42) 



and the inverse covariance matrix, 



/ 0.48435 -0.101383 
-0.101383 3.2882 
-0.164945 -2.45499 

-0.0305703 -0.0787898 

-0.097874 -0.252254 

V -0.106738 -0.2751 



-0.164945 -0.0305703 

-2.45497 -0.0787898 

9.55916 -0.128187 

-0.128187 2.78728 

-0.410404 -2.75632 

-0.447574 1.16437 



-0.097874 -0.106738 \ 
-0.252254 -0.2751 

-0.410404 -0.447574 
-2.75632 1.16437 

14.9245 -7.32441 
-7.32441 14.5022 / 



Finally we used the CMB shift parameter R 
#ov^mO Jo 1089 TTTFT- For tnis Xcmb is defined as, 



XCMB (#) 



(R(6) - i? ) 2 



(44) 



Where, Rq = 1.725 ±0.018 [36J. 

In this analysis we have two model parameters (3 and 
!l m o. We varied j3 from —0.5 to 1 and O m o from 0.2 to 
0.33. Figure (O shows the results in the (/?, Q m o) pa- 
rameter space. The analysis is performed for M = 10~ 3 
eV. Minimum value of the total \ 2 is at /? ^ 0.04 and 
^mo ~ 0.275. Using these best fit values of the model 
parameters growth index 7 is plotted in the figure © 
with the la and 2a errors. At late time, the evolution 
is consistent with Dark Energy models where w > — 1 
[37| . But the model shows a larger value of 7 for all 
redshift compared to ACDM. This is a different char- 
acteristic than /(i?)-gravity models [H| or scalar-tensor 
theories [39| where 7 < 0.55. 

We find that /3 ~ 0, /3 large is not a problem for local 
physics (chameleon, Vainshtein) but at large scales we 
need to consider j3 small as in coupled quintessence. In 
fact, the galilcon term is negligible at large scales and we 
recover the standard coupled quintessence model. 




-0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 



FIG. 5: The la (inside) and 2a (outside) likelihood contours 
in the (Q m , ,/3) plane for total XGro W th+SN+BAO+CMB- 



VIII. CONCLUSION 

In this paper, we have investigated a particular case of 
a more general class of models which do not produce a 
gravitational slip. This model includes the L 3 -term de- 
rived in the decoupling limit of DGP and a galileon field 
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ZBAO 


0.106 


0.2 


0.35 


0.44 


0.6 


0.73 




30.95 ± 1.46 


17.55 ±0.60 


10.11 ±0.37 


8.44 ± 0.67 


6.69 ± 0.33 


5.45 ±0.31 


D \( Z BA0) 



TABLE I: Values of „ , for different values of zbao- 




0.54 ; 

0.52 \ : 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

z 

FIG. 6: The la (inside) and 2a (outside) of the growth index. 
The central line is the growth index with the best fit values 
of the model parameters /3 and Q m o in the case of the double 
exponential (|10[) . 

potential added to the Lagrangian on phcnomcnological 
grounds. The potential breaks the galilcan symmetry 
in a flat spacctime but serves as an important toll at 
large scales for obtaining a viable cosmology. The scale 
in the model and the parameters are carefully set such 
that the galileon term is dominant at small scales which 
suppress the effect of the chameleon effect leaving the 
Vainshtein mechanism to operate. This was one of the 
motives of our proposal to replace the chameleon mech- 
anism as the latter is plagued by the problems of fine 
tuning and large quantum corrections due to the large 
mass of the chameleon field. In contrast, the domination 
of the galileon term at small scales successfully gener- 
ates a screening effect without requiring a large mass of 
the field and the quantum corrections remain small. We 
have, therefore shown that models including light mass 
scalar fields coupled to matter become viable as soon as 
we introduce a lower order galileon, term, (V^)) 2 D^> in 



the action. We should also emphasize that this term docs 
not significantly modify the background cosmology of the 
model. Our model is built in away as not to produce large 
deviation from standard model for cosmological pertur- 
bations due to the vanishing of the gravitational slip. 
This approach can be generalized to a large class of light 
mass scalar field models which are otherwise ruled out 
by local gravitational constraints. 

We have demonstrated that the model passes the local 
tests and can give rise to viable cosmology at late times. 
It produces a small deviation of the growth factor com- 
pared to ACDM model. It interesting to note that in case 
of pure coupled quintessence, G e // = G(l±2/3 2 ) thereby 
requiring coupling to be small in order to respect the 
BBN constraint. On the other hand, in the model under 
consideration, G c ff = G as in General Relativity which 
certainly alleviates the constraints on coupling at large 
redshifts. In principle, our scenario involves cosmological 
constraints at low redshift which depends on the form of 
the potential as usual. 
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